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ABSTRACT 

We performed, for the first time, the simulation of spiral-in of a star cluster formed close to 
the Galactic center (GC) using a fully self-consistent A^-body model. In our model, the central 
super-massive black hole (SMBH) is surrounded by stars and the star cluster. Not only are the 
orbits of stars and the cluster stars integrated self-consistently, but the stellar evolution, collisions 
and merging of the cluster stars are also included. We found that an intermediate-mass black 
hole (IMBH) is formed in the star cluster and stars escaped from the cluster are captured into a 
1:1 mean motion resonance with the IMBH. These "Trojan" stars are brought close to the SMBH 
by the IMBH, which spirals into the GC due to the dynamical friction. Our results show that, 
once the IMBH is formed, it brings the massive stars to the vicinity of the central SMBH even 
after the star cluster itself is disrupted. Stars carried by the IMBH form a disk similar to the 
observed disks and the core of the cluster including the IMBH has properties similar to those of 
IRS13E, which is a compact assembly of several young stars. 
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the SMBH is not clear. 



Young and massive stars have been found 
within one pa rsec fro m the Galactic cent er (GC) 



(|Krabbe et al . 1995; P aumard et"alll2006f) . Some 
of them are -1000 AU from the central SMBH. 
How these massive stars were brought to the vicin- 
ity of the SMBH has been a mystery. One possible 
scenario is the following. A star cluster formed at 
a few tens of pc from the GC, and then spiraled 
in due to the dynami cal friction j C erhard 2001). 
Previous simulation s ( Portegies Zwart et al]|2003l: 



tion in an accretion disk ( 


Levin & Beloborodov 


2003 


; Navakshin et al. 2007 


; Hobbs & Navakshin 


2008 


). Giant molecular clouds fall into the GC 



Kim fc Morrisl2003l : [Kkn et al.l2004HGurkan fc Rasi 



sic 



20051 ) have shown that the timescale of spiral-in 
of the star cluster can be short enough. How- 
ever, how close the stars can actually approach 



and form massive gaseous disks around the cen- 
tral BH. Stars form in the disk if it becomes 
gravitationally unstable and results in fragmen- 
tation. However, accretion disks have difficulty 
producing stars with eccentric orbits and a com- 
pact assembly of stars like IRS 13E, which is lo- 
cated at ~ 0.13 pc in projection from the GC 
and contains half a dozen young stars and prob- 
ably an IMBH dHansen fc Milosavlievidl J200I 



Maillard et all |2004 . iMapcUi ct al. (200^ 
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ar- 
gued that if a gas cloud undergoes a very close 
encounter (pericenter distance of 0.01 pc) with 
the central SMBH, the tidal compression could 
trigger the star formation, resulting in stars in 
close, bound orbits. However, how a cloud can 
come that close to the GC is not clear. On the 
other hand, if the star cluster has an eccentric 
orbit, the orbits of stars escaped from the cluster 
are also eccentric. The remnant of the core looks 
like IRS 13E. 
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We performed a fully self-consistent iV-body 
simulation in which the internal dynamics of the 
cluster, that of the parent galaxy, and interac- 
tions between cluster stars and galaxy stars are 
correctly handled. In previous simulations, when 
the internal dynamics was followed by an accurate 
A^-body code, the parent galaxy had to be mod- 
eled as a fixed potential with some fitting formulae 
for the dynamical friction. This means that the or - 
bital evolution is not accurate, ii et alj I (|2008l) 
showed that the actual orbital decay of the cluster 
is faster than that of previous simulations and the 
main reason is stars with a mass grater than 90% 
of an initial cluster escaping from this cluster. 

We describe the method of our A^-body simula- 
tion in section 2. In section 3 we show the results 
of simulations. Section 4 is for summary. 

2. Method 

2.1. Models and Initial Conditions 

We adopted two models for the galactic cen- 
ter. For both values, we adopted a model based 
on a King model with the non-dimensional central 
potential Wq = 10 as galaxy models. We placed 
the central SMBH with the mass of 3.6 x IO^A/q 
( Eisenhauer et al. l2005h and our galaxy models 
represent the central region of our Galaxy (see fig- 
ure [1]). The difference between two models corre- 
sponds to the initial position of a star cluster. One 
(GL) is for a run with the initial position of the 
cluster 12.5 pc from the GC, and the other (GS) 
5 pc from the GC. Their total masses on the real 
scale (excluding the SMBH) are 5.9 x lO'^M© and 
2.9 X IO^Mq, respectively. The number of par- 
ticles is 2 X 10^ for both models. GS has better 
mass-resolution, but a smaller half-mass radius of 
9.6 pc. GL has a larger half-mass radius of 22 pc 
for the farther initial position of the cluster. These 
models are summarized in table [U 

As a model of a star cluster, we adopted a King 
model with non-dimensional central potential of 
Wq — 6. For runs from the initial distance of 12.5 
pc from the GC, we used a model with 64k parti- 
cles and the total mass of 2.1 x IO^Mq (SC64k). 
For runs from 5pc, we used a model with 32k par- 
ticles and the total mass of 1.0 x lO^Mg (SC32k). 
Initial mass function of stars in the clusters is a 
Salpeter with lower an d upper cutoff at 1 and 
100 Mq (jSalpeteri Il955[ ). We assigned each star 



a mass randomly chosen from the Salpeter initial 
mass function, irrespective of its position. The 
tidal radii of the models are 1.1 pc for SC64k and 
0.65 pc for SC32k. They are smaller than the tidal 
limits at their initial positions. The cluster models 
are summarized in table [2l 

2.2. A^-body Simulation 



We used the Bridge code (jFuin et al.l 120071 ) to 
handle the interaction between the parent galaxy 
and the star cluster fully self-consistently. The 
Bridge scheme is a tree-direct hybrid scheme. 
Only the internal motion of the star cluster is 
calculated by the direct scheme with high accu- 
racy, and all other interactions are calculated by 
the tree algorithm. The splitting between the di- 
rect part and tree part is through the splitting 
of the Hamiltonian i n a way similar to the mixed 
variable symplectic (iKinoshita. Yoshida. fc Nakai 



19911 : IWisdom fc Holmanlll99ll ). With the Bridge 
scheme, we can treat a large-A^ system with em- 
bedded small-scale systems fully self-consistently 
and accurately. 

We used two sets of numerical parameters for 
each galaxy model. They are summarized in ta- 
ble [3l In our model, the softening length be- 
tween galaxy particles and the SMBH is 0.2 pc. 
Therefore, our simulation has the resolution limit 
around 0.2 pc for the motion of the star cluster 
and cluster stars within the parent galaxy. We 
used the opening angle 9 = 0.75 with the center- 
of-mass approximation for the tree. The Simula 



tion is performed using GRAPE6 (jMakino et al 
2003). 

2.3. Stellar Collisions and Evolutions 



In our simulation, we adopted collisions of stars 
in a star cluster and formation of an IMBH in 
the cluster. Recent simulations showed that in a 
dense star cluster, runaway collisio ns of stars form 

a very massive star CP ortcgics Z wart et al.l 1999t 

Portcgi es Zwart fc McMillanll2002i : |Portegies Zwart et al 
12004 : Jreitag et allbOOel T If it i s massive enough , 



it will collapse in to an IMBH (jFrver et all 12001 



Heger et al.l 120031 ). While the massive star grows 



through collisions, it loses its mass due to the 
stellar win d. Very massive stars lose their mass 



rapidly (jBelkus et al.l 120071 1. However, these pro 



cesses are not well understood for stars more mas- 
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sive than lOOOM©. Therefore, we treated the 
colhsion criterion and the mass-loss rate as pa- 
rameters. In this section, we describe our model 
of stellar collision and evolution and discuss the 
difference caused by these parameters. 

If the core collapse occurs in the star clus- 
ter, collisions between stars in the core occur 
rather frequently. We let two stars merge in 
the star cluster if they approach to twice the 
sum of stellar radii, which is also adopted in 
Portegies Zwart et al. ( 19991 ). If two stars ap- 
proach, the tidal capture occurs: the tidal force 
exerts stellar oscillations and the st ars lose their 
orbit al energy and becorne bina ry ("Fabian et al 



CAK model is given by 



ry I 

19751: iPress fc Teukolskvl 119771 : [McMillan et al 
19871 ). The critical distance for the tidal cap- 
ture is around 3-4 times the radius of the star 
for two identical stars. When a tidally captured 
binary is actually formed, it would merge fairly 
quickly because of perturbations from nearby 
st ars. For the radii of stars, we adopted the result 
of lHurlev et"aLl (|2000l ) (here after HPTOO). The 
structure of very massive sta rs with 100 < M < 
lOOOMp, is inves t igated by ilshii et al.l (|l999l ) and 
Yungelson et al. ( 20081 ). while that of stars more 
massive than IOOOMq is not. We used the results 
of HPTOO for very massive stars with > lOOM©. 
Mass ive stars (> lOOM fn) have a core-halo struc- 



ture (|lshii et al.l[l999l ). The result of HPTOO we 



adopt ed is close to the size of the core in llshii et al.l 
(Il999l) . Therefore, we adopted twice the sum of 
the radii as the collision criterion. To see the effect 
of the different criterion, we also performed sim- 
ulations using the s um of the stellar radii as the 



collis i on criterion (IPortegies Zwart fc McMillan 
' iFreitag et al 



Portegies Zwart et 



all 120041: 



2002; 
20061) . 

We also take into account the mass loss 
due to the s t ellar wind for very massive stars. 
Belkus et al.l (|2007l ) have investigated the mass- 
loss rate for massive stars with < IOOOM0. How- 
ever, the stellar evolution of stars with > lOOOM© 
is not well understood. Furthermore, merged stars 
show a evolut ion different from that of single stars 
(jSuzuki et al.li2007) . Therefore, we treat the mass- 
loss rate as a parameter. 

Our model for the stellar wind is based on 
line-driven winds developed by Castor et al. 
(1 975; hereafter CAK ). We used the formalism 
of lOwocki et"all tOOi ). The mass-loss rate of the 



1 — ct 



QTe 

1-r, 



(l-a)/a 



(1) 



where is the luminosity of the star, Fg is the 
Eddington parameter, and a and Q are the power 
index and normalization of the line opacity distri- 
bution. The Eddington parameter is given by 



inGA'L 



(2) 



where M, is the mass of the star, G is the grav- 
itational constant, c is the speed of light, and 
is the opacity. For fully ionized plasmas with 
hydrogen mass fraction X, the opacity is given 
by K = 0.2(1 + X) cm^ g~^, where we assumed 
X = 0.7. For very massive stars, the luminos- 
ity is proportional to the mass because in such 
massive stars, the contributi on of radiation to th e 
total pressure is very large ( Marigo et al. 20031 ). 
Hence, we assumed M* oc M*, here M* is the 
mas s of the star. W e adopted a = 0.5(CAK), Q — 
10^ (|rxavlevlll995l ). = 3.2 x 10^{M^/Mq){Lq) 
(Suzuki, private communication) , which is close to 
the Eddington luminosity. From these values, we 
obtained 

Af, = 9.66 X 10-« (^^) (Mo/yr). (3) 



On the other hand, the recent results (Belkus et al 



20071 ) show a higher mass-loss rate. Therefore, we 



also adopted a mass-loss rat e five t imes higher, 
similar to the result of Belku s~et al.l (I2OO7). The 
evolution of the mass for stars with 1000 Mq ini- 
tial mass is shown in figure |2l We neglected the 
mass loss of the stars less than 300 M© because 
they are small enough. 

We assumed that very massive stars formed by 
runaway collisions form black holes ( BHs) at the 
end o f their main - seque nce lifetime ( Frver et al.l 
2OOII: iHeger et all 12003^ For the lifetime, we 



adopted the results of HPTOO for stars with masses 
less than 100 Mq, interpolated the results of 



Belkus et al.l (|2007l ) for stars with masses up to 



1000 Mq, and extrapolated them for more than 
1000 Mq. We assigned new ages to stars born by 
stellar collisions using th e follo wing formalism of 
iMeurs fc van den Heuvell (|l989l ): 



<ago(™l+™2) = : ^ iago(™l)' \V 

mi + m2 T^s(mi) 
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where toi and rn2 are the masses of stars (mi > 
W2), iago is the age of the star, and Tms is the main- 
sequence hfetime. The new star looks rejuvenated. 
However, the assumption used to derive equation 
Q is that the convection core of the primary will 
occupy the core of the merger with the mass and 
com position unchanged. [ Suzuki et alJ pOOTH and 
also Icaburov et ah ( boosh found that in the case 
of unequal-mass merger the core of the secondary 
will first occupy the center of the merger, resulting 
in significant mixing-in of hydrogen into the new 
convection core. So equation (4) might underesti- 
mate the lifetime of the merger. We assumed that 
the massive star directly collapses to an BH. The 
final evolution of the very massive star depends on 
its helium core mass at the end of its life. If it is 
more massive than ~ ISOMf:^ it collapses direc tly 
to an BH without mass loss ( Frver et al. 200l[ ). 

We performed simulation for isolated star clus- 
ters and investigated the difference caused by the 
parameters. The model of star cluster is SC64k. 
The runs are summarized in tablelH Each parame- 
ter has two values. The results are shown in figure 
[3l The difference due to rcou is smaller than that 
due to the mass-loss rate because the stellar radii 
is smaller than the radius, but it becomes larger 
when the mass grows. The mass-loss rate signifi- 
cantly affects the final mass. In any case, however, 
the mass accretion rate due to stellar collisions was 
higher or almost the same as the mass-loss rate. 
The mass loss will not prevent the formation of an 
IMBH in the star cluster. 

2.4. Runs 

We performed fully self-consistent iV-body sim- 
ulations of a star cluster and its parent galaxy sys- 
tem. Runs are summarized in table El For the 
mass-loss rate and collision criterion, we adopted 
two extreme cases. They are upper and lower lim- 
its of the mass of IMBHs. 

3. Results 

3.1. Evolution of Star Cluster and IMBH 

The orbit of the star cluster decayed due to the 
dynamical friction and the cluster was disrupted 
by the tidal force. The evolutions of star clusters 
are shown in figure |4l Top and middle panels show 
the orbital and bound-mass evolution of the star 
clusters. Bottom panels show the mass of the most 



massive star in the star cluster. Black points in 
the bottom panels show the times when IMBHs 
formed. 

In each case, when the pericenter distance be- 
comes less than 2 pc, the star clusters are almost 
completely disrupted. The disruption time, rdis, is 
summarized in table [6l We defined the disruption 
as the time when the number of stars bound to 
the cluster becomes less than ten. 

Before the cluster was completely disrupted, an 
IMBH formed in the cluster through the runaway 
collisions of stars. A massive star grew via re- 
peated collisions and finally the mass of the mas- 
sive star reached 3000-16000 Mq. The evolution 
of the massive star is shown in bottom panels of 
figure m The IMBH mass strongly depends on the 
mass-loss rate. The time when the IMBH formed 
was earlier than the disruption time except for the 
run HS32k. These results are summarized in table 

151 

After the disruption of the star cluster, the 
IMBH formed sinks into the GC on the time 
scale of the Chandra sekhar dynamical friction 
(jChandrasekhaii Il943l ). The time scale of the 
orbital evolution due to the dynamical friction 
is proportional to the mass of the star cluster 
(IMBH). The IMBH mass of LD64k is around 
three times more massive than that of HS64k. 
The orbital evolution of LD64k is roughly three 
times faster than that of HS64k. On the other 
hand, the orbital decay of the cluster looks faster 
than the time scale. This is because the dynami- 
cal friction was e nhanced due to th e stars escaping 
from the cluster ( Funi et al. 20061 ). 

The orbital evolution of the star cluster became 
slower and slower at the end of the simulation. 
This is because we adopted a large softening length 
between the SMBH and galaxy particle, which is 
0.2 pc. If the softening length were smaller, the 
cluster would approach closer to the GC and carry 
resonant stars. We will perform a simulation with 
a smaller softening length and report the results 
in a forthcoming paper. 

Figure [5] shows the projected distribution of 
stars at the end of the run (7.25Myrs) of LD64k. 
Only stars which were originally members of the 
star cluster are plotted. The arrows in the right 
panel show the proper motions of the stars. Red 
ones move clockwise and green ones move coun- 
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terclockwise. This figure lo oks very similar to K- 
band images of the GC (cf. ILu et ai. 2006L figure 
2). In particular, they reached a distance of ^ 0.2 
pc from the GC. We can see that, even though 
the cluster is almost completely disrupted at 2pc, 
there are a number of stars which were brought 
much closer to the central SMBH. 

3.2. Trojan Stars 

We found that many stars were carried near 
the GC though almost all stars became unbound 
at around 1-2 pc from the GC. Figure [6] shows how 
these stars were brought near the central SMBH. 
Red solid curves show the orbit of the star clus- 
ter. Here, we plotted the trajectory of six stars of 
LD64k. All of them have escaped from the clus- 
ter by 3Myrs, as can be seen from the separation 
of black and red curves. One would imagine that 
their orbits do not evolve once they have escaped 
from the cluster. However, actually these stars 
"follow" the spiral-in of the cluster remnant. This 
behavior occurs due to the 1:1 resonance with the 
IMBH formed in the cluster. These stars have a 
pericenter distance of less than O.lpc at the end of 
the simulation. 

In figure [6l stars R1-R4 were brought to the 
GC by the 1:1 mean motion resonance. Figure [7] 
shows the orbit of Rl in a rotating frame where the 
distance between the SMBH and IMBH is scaled 
to unity. Star Rl escaped from the cluster at T ~ 
2.5 Myr and was orbiting around the SMBH until 
T ~ 3.8 Myr. At T ~ 3.8 Myr, it was captured 
into the resonance. Its orbit was a horseshoe orbit 
for T ~ 3.8 Myr to T ~ 5 Myr. From T ~ 5 
Myr to T ~ 6.5 Myr, the orbit of Rl is that of 
a retrograde quasi-satellite. After that, it escaped 
from the resonance. Stars R2 through R4 show 
similar behaviors. Their distances from the GC 
have become smaller by a factor of five or more 
while they are in the resonance with the IMBH. 
Some of them escape from the resonance on the 
way to the GC. 

Some stars are scattered by the IMBH during 
the resonance and change their orbit into those 
approaching the GC. We categorized such stars 
into a Resonance and Scatter group (RS), and 
show one example in figure [6l In this case, its 
semi-major axis docs not decrease gradually, but 
suddenly through close encounters with the IMBH 
(see figure m middle-bottom panel). These reso- 



nant stars (R and RS) make up more than 90 % of 
stars whose pericenter distances, Tp, are less than 
0.5 pc from the GC (Table El). 

The right-bottom panel in figure [5] shows the 
orbital evolution of a star which escaped from the 
star cluster by a slingshot occurring in the star 
cluster. We categorized such stars into a Sling- 
shot group (S). Some of them come very close to 
the GC. Moreover, these stars have very high in- 
clinations (some have retrograde orbits) to the or- 
bital plane of the star cluster because the escape 
direction from the star cluster is isotropic. They 
may be the origin of the stars which have high 
inclinations in the GC. Table [7] summarizes these 
results. We can see that the 1:1 resonance with 
the IMBH is the main mechanism which brings 
the stars close to the GC. 

Figure [8] shows the cumulative number of es- 
caped stars as a function of pericenter distance 
at the end of simulations. Solid curves include 
all stars and dashed curves only massive stars 
(> 20Mq). Dotted curves show the expected num- 
ber of stars if the fraction of stars is independent of 
Tp . We can see that there is a much larger number 
of massive stars in the central sub-parsec region 
than expected from the initial mass function. 

In the case of run LD64k, there are 16 mas- 
sive stars (m > 2OM0) with the pericenter dis- 
tance less than 0.5 pc. Five of the 16 stars reached 
within 0.2 pc and three reached within 0.1 pc (see 
table [H). More than 10 % of stars with Vp < 0.5 
pc have masses 2OM0. If the distribution of stars 
does not depend on their masses, the expected 
number of massive stars for the radius of O.lpc 
is around 0.3. If there is no mass-dependent evo- 
lution effect, the chance that we find more than 
three stars with rp < 0.1 pc is lO^'^. In the other 
cases, a few persent of stars near the GC was mas- 
sive at Tp = 0.5 (pc) and the fraction is higher 
at more deeper region. In the initial mass func- 
tion, only 0.5% stars have mass more than 20Mq 
and ~ 0.1% after runaway collisions because many 
massive stars are used up for the formation of the 
IMBH. 

Figure [8] also shows that the number of the Tro- 
jan stars depends both on the mass of IMBH and 
the spiral-in timescale of IMBH. The number of 
resonant stars decreases when the IMBH mass is 
smaller. However, as the IMBH spirals in towards 
the GC in 5-6 Myrs, the number of Trojan stars is 
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not much different from tliat in our standard run 
(LD64k). 

These results imply that star clusters selectively 
carry massive stars close to the GC. Massive stars 
sink to the center of their parent star clusters due 
to the mass segregation, while the stripping by the 
tidal force removes less massive stars from the out- 
skirts of star clusters. As a result, massive stars 
tend to remain in star clusters and are carried to 
a few parsec from the GC before they escape from 
the cluster. Slingshot stars also include many mas- 
sive stars because slingshots occur at the center of 
the star cluster where massive stars are gathering. 
Thus, star cluster scenarios can naturally explain 
why massive stars lie within the central parsec. 

3.3. The Remnant of the Core and IRS 
13E 

At the end of the simulation, only one mas- 
sive star is bound to the IMBH. However, sev- 
eral stars were bound till T = 6 Myrs for LD64k. 
These bound stars look like IRS 13E, which is lo- 
cated at ~ 0.13 pc in projection from the GC 
and contains half a dozen massive and young stars 
within ~ 0.01 pc. Because of their very similar 
proper motions, it has been suggested that they 
are bound to an IMBH and IRS 13E is the rem- 
nant of a star cluster containing an IMBH with 
~ 10^ - lO^ Mrr. dHansen fc MilosavlievicI 120031: 



Maillard et al.l 120041 ) . Figure [9] shows the proper 
motion of stars bound to the IMBH at T = 4.62 
Myr for LD64k. The tidal radius is 0.15 pc and 
larger than the frame, but a few stars are within 
0.01 pc. This is similar to the observation o f the 
IRS 13E (see figure 2 in lMaillard et al.ll2004f ). 



4. Summary 

Using A^-body simulation, we showed that 
many young and massive stars are carried to the 
GC by a star cluster due to the 1:1 mean motion 
resonance with an IMBH which is formed in the 
cluster. In addition, we found that slingshots in 
the star cluster throw stars into orbits which pass 
near the GC. These orbits have very high inclina- 
tions and are sometimes retrograde orbits. They 
are new channels which carry young stars to the 
central parsec. Our simulation demonstrated the 
existence of massive stars and we explained why 
they form a disk-like structure. The possible exis- 



tence of two counter-rotating disks might suggest 
that two clusters have spiraled in. The interac- 
tion and resonance between stars and IMBHs from 
multiple clusters might be responsible for the ex- 
istence of stars which are very close to the central 
SMBH. 
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Table 1: Models for the galaxy 





King Wo 


N 


MG(Me) 


m(MQ) 


Mbh(M0) 


^■h (pc) 


GL 


10 


2 X 10" 


5.8 X 10^ 


29 


3.6 X 10" 


22 


GS 


10 


2 X 10^ 


2.9 X lO'' 


15 


3.6 X 10^ 


9.6 



Table 2: Models ior ihc star cluMcr 



King Wo 




Msg (Me) 


Tc (pc) 


Th (pc) 


n (pc) 


SC64k 6 


65536 


2.0 X 10° 


5.9 X 10"^ 


1.6 X 10-1 


1.1 


SC32k 6 


32768 


1.0 X 10^ 


3.5 X 10-2 


9.6 X 10-2 


6.5 X 10-1 



Table 3: Parameters 





esc (pc) 


ec (pc) 


esc-BH (pc) 


ec-BH (pc) 


Aitree yr 


PI 


0.0 


4.9 X 10-^ 


2.5 X 10-^ 


2.5 X 10-1 


2.4 X 10^ 


P2 


0.0 


2.3 X 10-2 


2.4 X 10-2 


2.4 X 10-1 


1.1 X 102 
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Table 4: Runs for isolated star clusters 





m 


»'coU 


MiMBH(Af0) 


Mlost(MQ) 


LD 


Low^ 


Double=^ 


1.8 X 10-* 


2.9 X 10^ 


LS 


Low 


Single^ 


1.4 X 10^ 


2.0 X 10^ 


HD 


High^ 


Double 


8.3 X 10-'' 


9.0 X 10^ 


HS 


High 


Single 


6.7 X 10^ 


6.8 X 10^ 



im = 9.7 X (Mo/yr) 
= 5.0 X lO-'^m (Mq /yr) 
^'■coll = ri + r-z 
^^coll = 2(ri +r2) 

Table 5: Runs 





Galaxy 


Star cluster 


m 


T'coU 


Parameter 


Ro (pc) 


LD64k 


GL 


SC64k 


Low^ 


Double^ 


PI 


12.5 


HS64k 


GL 


SC64k 


High^ 


Single* 


PI 


12.5 


LD32k 


GS 


SC32k 


Low 


Double 


P2 


5 


HS32k 


GS 


SC32k 


High 


Single 


P2 


5 



= 9.7 X IQ-^m {Mq/jx) 
= 5.0 X IQ-'^m (Mq/yt) 

^rcoll = ri + r-2 

^^coll = 2(ri + r2) 

Table 6: Results 





iiMBH(Myr) 


Mimbh(M0) 


Afiost(Af0) 


tdis' (Myr) 


LD64k 


2.6 


1.6 X 10* 


2.4 X 10-^ 


5.3 


HS64k 


2.6 


6.3 X 10^ 


6.5 X 10^ 


4.3 


LD32k 


2.3 


8.8 X 10^ 


1.5 X 10^ 


2.7 


HS32k 


2.2 


3.0 X 10^ 


4.1 X 10^ 


1.7 



^Disruption time of star clusters, when the number of bound stars becomes less than ten. 
Table 7: The number of stars within 0.5 pc for LD64k. 



Pericenter distance 


Bound 


Slingshot Resonance 


Resonance & Scatter 


total 


Tp < 0.5 pc 


2 


27 280 


45 


354 


rp < 0.2 pc 





4 52 





56 


Tp < 0.1 pc 





23 





23 



Table 8: The mmibcr of massive stars within 0.5 pc for LD64k. 



Pcricciiti^r (listaucc 


Bound 


Sliiigsliot, 


R(>sonancc 


Rt^sonancc & Sea 


.tor 


Total 


Tp < 0.5 pc 


1 


5 


7 


3 




16 


r,, < 0.2 pc 





1 


4 







5 


Tp < 0.1 pc 








3 







3 
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R (pc) 



Fig. 1. — Density profiles for our models and the Galactic center (|Genzel et al.ll2003f) . 




r(Myr) 

Fig. 2. — Evolution of the mass for stars with initial mass = lOOOM©. The solid and dashed curves are the 
results with low mass-loss rate (see equation Q in the text) and high mass-loss rate (five times as that of 
equation ([3])), respectively. 
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12 3 

T (Myr) 

Fig. 3. — Evolution of the most massive star in the star cluster. The model of star cluster is SC64k. Black 
points show the time when IMBHs formed. "L" and "H" mean low- and high-mass loss rates, and "D" and 
"S" denote the coUision radius. See tabled] for details. 
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T (Myr) T (Myr) 

Fig. 4. — Evolution of star clusters. Top and middle panels show the orbital and bound-mass evolution of 
star clusters, respectively. Bottom panel shows the evolution of the most massive star in the star clusters. 
The black dot shows the time when the massive star became an IMBH. 
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Fig. 5. — Projected distribution of stars at T = 7.26 Myr for LD64k. The orbital plane of the star cluster 
has an inclination oi i ~ 127° with respect to the plane of the sky, w ith a half-line o f ascen ding nodes at 
n = 99° east of north. These values are used to mimic the result of Paumard et al. I llooi). Right panel 



shows the central region shown by a red square in left panel. Arrows in right panel show the proper motion 
of stars. Red one shows clockwise orbit and green shows counter-clockwise. Red and white crosses show the 
positions of the IMBH and SMBH, respectively. 
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Fig. 6. — Evolution of position, eccentricities, and inclinations of stars escaped from the star cluster for 
LD64k. Red curve shows the orbital evolution of the star cluster or the IMBH formed in the cluster. Green 
and yellow lines show the time when the IMBH was formed and when the star cluster was disrupted (less 
than ten stars are bound to the cluster), respectively. 
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X X 

Fig. 7 — Orbit of star Rl in a rotational frame, where the IMBH is fixed at (1.0, 0.0). The SMBH Hes on 
the origin. The distance of the star from the SMBH is normaUzed by the distance between the SMBH and 
the IMBH. Gray dots show all positions of the star obtained from snapshots. Colored curves or dots show 
the orbit or position of the star within time spans given in each panel. 
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rp (pc) (pc) 

Fig. 8. — Cumulative number of escaped stars as a function of peri-center distance at the end of simulations. 
Solid curves include all stars and dashed curves only massive stars (> 20Mq). Dotted curves show the 
expected number of stars if the fraction of stars is the same, irrespective of Vp. 
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Fig. 9. — Snapshot of the remnant of the core. V ectors show the prop er motion of stars for massive stars 
(> 2OM0). It is similar to IRS 13E (see figure 2 in iMaillard at al.ll2004l ). 
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